Detection of Spin Correlations in Optical Lattices by Light Scattering 



Ines de Vega, 1 J. Ignacio Cirac, 1 and D. Porras 1 

1 Max- Planck- Instttut fur Quantenoptik, Hans-Kopfermann-Str. 1, Garching, D-85748, Germany. 

(Dated: February 2, 2008) 

We show that spin correlations of atoms in an optical lattice can be reconstructed by coupling the 
system to the light, and by measuring correlations between the emitted photons. This principle is 
the basis for a method to characterize states in quantum computation and simulation with optical 
lattices. As examples, we study the detection of spin correlations in a quantum magnetic phase, 
and the characterization of cluster states. 
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Ultracold atoms in optical lattices open exciting 
prospects for the investigation of quantum many-body 
phases in a highly controllable setup. For example, spin- 
spin interactions between atoms in a Mott insulator state 
can be tuned in such a way that one can simulate a rich 
variety of models from quantum magnetism [l[ . Further- 
more, this system is ideally suited to implement a quan- 
tum register for quantum computation ,2] , and to realize 
multiparticle entangled states with cold controlled colli- 
sions 0, An important benchmark in this direction 
has been the creation of a cluster state since it can 
be used as a resource for measurement based universal 
quantum computation 0, 0] • 

A major disadvantage of this setup is the fact that 
atoms are separated by optical wavelengths, and thus it is 
difficult to address them individually. This impedes one 
to measure directly the spatial dependence of spin-spin 
correlations, which is essential in the characterization of 
quantum phases. For these reasons the development of 
accurate methods to measure properties of atomic op- 
erators is basic to the usefulness of optical lattices for 
quantum computation and simulation. A possible solu- 
tion relies on the detection of the atom number distribu- 
tion in time of flight (TOF) imaging. Within this tech- 
nique atomic density operators in momentum space are 
measured by taking absorption images of the expanding 
cloud, after having switched off the trapping potentials. 
From TOF images density-density correlation functions 
in momentum space [H, [3, [Toj ] are reconstructed by de- 
tecting statistical correlations between different images. 
On the other hand, one may use quantum non-demolition 
techniques on the atom-light interaction in order to de- 
termine certain atomic collective observables 
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In this letter we describe an alternative approach which 
allows one to measure spin correlations without switching 
off the optical potential. Furthermore, it is not necessary 
to drive the atoms out of the strongly correlated regime, 
thus avoiding certain complications that are met in TOF 
measurements of interacting fermions. Our method is 
based on off-resonant scattering of an incident laser with 
the trapped atoms. Scattered photons with a certain po- 
larization are coupled to certain spin operators. In that 
way correlations between photons emitted in different di- 
rections are proportional to ground state correlations of 



spin operators in momentum space. Different detection 
schemes, such as photon counting or homodyne detec- 
tion, may give information of different types of correla- 
tions. One of the strengths of the method is that it is 
based on photon detection, and not on atom detection. 
For this reason, detection of any possible spin-spin corre- 
lation is achieved here naturally, just by controlling the 
polarization of the lasers and of the detected photons. 
The correlation of an arbitrary number of spin operators 
can be achieved by considering correlations of different 
detections. 

We consider an optical lattice filled with N atoms 
in a Mott insulator state with at most one atom per 
site. Each atom carries an arbitrary ground state hy- 
perfine spin J, and has an optical dipole transition of 
frequency u>q to an excited states manifold of spin J'. 
This transition is coupled with an off-resonant laser of 
frequency ojl, wavevector k^, and a certain linear polar- 
ization (7£. Photons are then scattered with momentum 
k and two different polarizations j3, such that (3 ^= o^. 
Light modes with momentum k are coupled to atomic op- 
erators J^ k = l/v^Ej jf,e iAk - r 3, with Ak = k - k L , 
where Jg are atomic spin operators, and the vector de- 
notes the position of the atoms in the lattice. This allows 
us to describe correlations of spin operators with differ- 
ent momentum, which as it will be later discussed pro- 
vides us in principle with full information to reconstruct 
spin correlations in the position space. We present the 
following results: (i) By measuring correlations between 
photons with momenta k we reconstruct correlations be- 
tween atomic spin operators J^ k , such that k = kh- 
Full access to correlations both in momentum and po- 
sition space requires that /cl^o > V^k (VStt) in 2D (3D) 
square lattices, where ro is the lattice constant. How- 
ever, these conditions are not always necessary to char- 
acterize a state in the momentum space, (ii) Measured 
correlations are in the ground state provided that the 
measurement time T -C where T is the spontaneous 
emission rate. The validity of this approximation is stud- 
ied by simulating the emission properties of interacting 
atoms in a lattice by means of a bosonic description of 
spin excitations, (iii) As an illustration of applications, 
we discuss the measurement of a magnetic phase, an ex- 
ample that is relevant to quantum simulation, as well as 
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a measurement of a cluster state in a 2D lattice, which 
has an important application in quantum computation. 

The incident laser has a large detuning A = ujj/ — u>l, 
with respect to the uppermost level of J', with frequency 
ujjr. In that situation, the excited manifold can be adi- 
abatically eliminated [LJ, [l4|, and the evolution of any 
operator A acting on the spin J manifold, can be written 
in terms of an effective Hamiltonian for light -matter in- 
teraction (d/dt)A — i/H[H^, A}. Along the same lines as 



Heisenberg equations of the field operators, we find 
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where the sum goes over a',/3' = x,y,z, and we have 
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|, the effective Hamiltonian can be written aennea J\n , ~ l^p e l^a',a k/} \, wnere e is 
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where g = d 2 ai/A, d is the atomic dipole matrix 
element, and a\ is a constant that depends on the 
particular transition considered. The electromagnetic 
field is usually expressed in terms of the two or- 
thonormal polarization vectors e k ;_l_k (/ = 1,2) and 
the corresponding creation (annihilation) operators 
(aki)- Nevertheless, it is more convenient for us to 
express Ej ot in the laboratory frame. We decompose 



the field in the following way: E^ ot = JZ^E^" 
W a + +E{ + ), where E£+ = £ k e k a k(T e- lA ^ +lk ^<7, and 
E^+ = J2k e k a k<r e ~ lAfc t+lk r 3(T are fast and slowly rotat- 
ing terms respectively 15 1. and E^ + = e£e -lkL ' Tj &l is 



the laser field. The constants are e k = WhuJk/ (2eou), 
= cjfe ± Wi, and a^, a = x,y, z. In addition, we have 
defined photon operators in the laboratory frame, a k:E — 
cos 9 cos </>a k i + sin 0a k 2 , a kj , = cos 6 sin </>a k i — cos 0a k 2 
and a kz = — sin #a k i , with 6 and the angular coordi- 
nates of the wave vector k in the laboratory frame. 

Let us consider as an example the H eS when the laser 
polarization is ul — z, 



H° s = igJ2[Ji(E; : 



-) - 4(Et - E?)] (2) 



where E 3 ~ = E 3 ' e~ 1 ^ ^ - E%~ e***-** , for a = x,y. 
The Hamiltonian © describes a coupling between the 
emitted j/-polarized (x-polariced) photons and the spin 
operators J 3 X (Jy)- Different laser polarizations give rise 
to the scattering of photons that are coupled to other 
spin operators. 

We show first how the light-matter coupling ([2]) allows 
us to measure equal-time spin correlations of the kind 
(J A %- Ak >, which are very useful in the characteriza- 
tion of many-body spin phases. As we show below they 
are related to N^(T) = v/(2ir) 3 J dfcfc 2 (a kct (T)a k/3 (T)), 
where v is the quantization volume. The diagonal ele- 
ments of this quantity are the number of photons emitted 
in the direction k during a time T, whereas nondiagonal 
terms may be readily obtained by rotating the polariza- 
tions of the photons prior to the measurement. From the 



the Levi-Civita symbol. A few approximations can be 
considered in order to simplify the expression ©. First, 
the evolution of the atomic correlation due to the atomic 
Hamiltonian i? at occurs in a time scale Ta = 1/e (where e 
is a typical eigenenergy of H at ), that is much larger than 
the light-matter interaction processes that are here con- 
sidered, and can therefore be neglected. Second, T can 
be made short enough as to ensure that the dependency 
of the correlation over r can be neglected. Finally, the 
condition T ^> tq (where tq is the environmental decay- 
ing time), allows us to extend to infinity the integration 
limits of the first integral. Hence, we can write 
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where the constant To = g 2 (^r) ttcl J dkk 2 e 2 .5(k — hi) 
is the spontaneous emission rate of a single atom. More- 
over, if T <C 1/r, where 1/r is an estimate of the spon- 
taneous emission decaying time of the system that may 
be renormalized by collective effects, we can write 



K p {T) « NTo ^ M* a ,M k m T(J**(0) J^(0)). (5) 

a' 13' 

The validity of this approximation will be later studied 
in more detail, since the measuring time T has to be 
also long enough as to ensure that a sufficient number of 
photons to characterize the state is detected. 

Apart from the photon number, other operators like 
the field quadratures are linked to spin observables. Let 
us consider detection in the far field limit, so that the 
detector is placed at a position R such that k^R 1- 
In that case, it can be shown that the modes that con- 
tribute more to the field are those with wavevector in 
directions k = R, with R a unit vector in the direction 



of the detector [17| . In addition, due to energy conserva- 
tion, the largest contribution to the emitted field comes 
from k — ki. Hence, the positive component of the emit- 
ted field corresponding to the a polarization ca be writ- 
ten as E+(r,T) « VN^e-^ t+ ^ R J2 p Mf a J^ AR (0), 

where AR = &lR — k^. Here we have followed simi- 
lar approximations as in the derivation of ((4J, and also 
that J ( Ak (r) « J^ k (0) for re[0,T], provided again that 
T <C 1/r. By performing homodyne detection one may, 
for example, measure the quadrature X k = e k (a kQ +a kct ), 
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which is related to the spin-operators in the following 
way: 

X*(T) = 2VA^]T.M k /3 J Ak (0), (6) 

9 p 



where k = fc^R, and j3 = x,y, z. We have de- 
fined J Ak (0) = (J Ak (0) + JJ Ak (0))/2 = £\ cos(Ak • 
Tj)Jp/y/N. From ([6j> , we see that a certain quadrature 
is proportional to a combination of system spin opera- 
tors. If we want to detect only a certain spin component, 
for instance J Ak , we should fix the detector in 9 and <j> 
such that only A4 ax ^ 0. Different values of Ak are 
scanned by changing k^, what would provide us with in- 
formation of atomic spins (and their correlations) within 
a 2D slice in momentum space. Spin operators in the 
whole 3D momentum space, may be obtained by con- 
sidering measurements in which both the detector and 
the laser are moved. Note that to extract information 
on spin operators from optical measurements requires to 
invert the matrix .M k = {A^ k g}, what in general can be 
done except for some particular values of 9 and 4> such 
that det[M] ^ 0. 

We now turn to study further the approximation 
N a ,f3{T) ~ N at p(T), which allows us to reconstruct 
ground state correlations from the emitted photons. The 
idea is to calculate how much of the information we get 
from the atoms corresponds to their ground state. The 
study of the radiative emission of an ensemble of inter- 
acting atoms poses a complicated many-body problem. 
We address it here by considering that atoms have a well 
defined magnetic ordering along the z axis, such that 
[J|, Jy] — iSjj, in a Holstein-Primakoff (HP) approxi- 
mation. Our detection method relies on being small 
up to T such that a sufficient number of photons has been 
emitted to ensure a good detection efficiency. Within the 
HP approximation, Hamiltonian ([2]) yields the following 
system of evolution equations: 



dt 
dt 



= -uj rffi k ^/- Ak [M^iJf) + M k x (J Ak >] , 
= uj dn k s«+ L Ak [M v k „(J- Ak ) + M y k E (J- Ak f ) 



where U — Ag 2 To, and the spin dynamics induced by 
iJ a t during the radiative time scale has been neglected. 
We have also assumed that atoms are distributed in a 
lattice with single occupation, something that does not 
alter the qualitative features of the emission process. 



'l/L 



J2j e l(1 ' rj is a function centered at q = with 
width 1/L. In the limit 1/L <C (1q, spin operators vary 
much more slowly than S^, L in momentum space, such 



r ±Ak\ 



that one can approximate J dQkJ\4^6^f L k (Lg / ~ 

(jj q ) / dtluM*^^, and Q becomes a closed sys- 
tem for each atomic operator with momentum q. Finally, 



from the quantum regression theorem, ((7|) can be used to 
evolve two operator averages, and calculate the emission 
pattern by means of Eq. (fj|. 

To present a definite situation, we consider atoms with 
spin 1/2, with H at = BJ2j J i + J J2<ij> J l J i' where 
B is an external magnetic field, we consider a ferro- 
magnetic interaction (J < 0), and < ij > refers to 
nearest neighbors. In the regime B > |J|, and far 
from the critical point, the system is in a paramag- 
netic phase with spins aligned along the z axis. The 
HP approximation yields an adequate description of the 
quantum fluctuations in the ground state. Within this 
approximation the problem can be rewritten in mo- 
mentum space, H at = £ q (2L? + J q )J^J~ q + JgJy' 1 , 
where we have used the Fourier transformed spin-spin 
interaction, J q = 2J(cos(g x r ) + cos(gj,r ) -I- cos(q z r )). 
Note that the description of the many-body problem in 
terms of spin- waves suits perfectly our purpose, since 
it also allows us to describe the emission process. The 
ground state correlations within this approximation are 



J«Jp) = ^T+JJIB, (jqj-i) = 1/2^1 + J q /4£, 



i/2, and (J q J~ q ) = -i/2. They are used 
as initial condition to the set of equations obtained by 
means of the evolution equations [[?]), 

As an example of measurement, we show in Fig. (TT]) 
the quantity Af^y'^ , for a certain small T such that 
most of the photons come from the system ground state. 
In addition, the relative error, 

(L^) 2 = J dco S (9)(<^-N^) 2 / J dco S (9)(Af^) 2 , 

(8) 

where we have fixed <fi = 0, is plotted in Fig. (fTJ) with 
respect to the number of emitted photons (proportional 
to T). It can be seen that the error remains relatively 
small when enough number of photons (of the order of 
600) have been emitted. We stress that the radiative 
emission from an interacting spin system is indeed an 
interesting problem by its own, which shows the interplay 
between quantum dissipation and many-body effects. 

The information gathered by measuring in different 
directions, may be used to reconstruct atomic correla- 
tions in the position space. To study the conditions 
which are required for this task, we focus on the par- 
ticular example of a 2D optical lattice. This case is rele- 
vant to the characterization of cluster states and exper- 
iments which simulate high-Tc superconductors. Con- 
sider also for concreteness that the incident laser propa- 
gates along the z direction, and that the lattice is de- 
fined within the x-y plane. Following ([6]), an homo- 
dyne detection of the quadrature X£ allows one to mea- 
sure the operator J Ak , where now we have that Ak = 
/cl (sin 9 cos 4>, sin 9 sin</>, cos 9 — 1). Then the spin oper- 
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(n,m) T (n,m) T (— n, — m) . ... 

- J R + J)j m the position space, 



'cp - ' J ^ ' J fi _ _ 
with (n, m) denoting a particular site within the lattice, 
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FIG. 1: Relative error between Nyy'^° and N yy e '* =0 (T), 
with respect to the number of emitted photons with a — y, 
integrating all over 6 for <f> — 0. Triangles: J/B = —0.5, 
Boxes: J/B = -0.1, Stars: J/B = -0.01.; Inset: Number 
of y-polarized photons for each value of 9 with = 0, with 
T = 0.001 (units of T = 1), and J/B = -0.5. Triangles: 
Myy^ = F N r o T dt{jf'^°(t)jf^=°(t)) (Emitted pho- 
tons). Squares: iv" y fl -*= (T) = r o iVT(j|^ =o (0) J*'*=°(0)) 
(Emitted photons that correspond to ground state) 

can be obtained as an inverse Fourier transform of J A £, 

r-l _ pi 

j'o'p = C J d~k x J dky cos(k x nn + k y m7r) Jq^ , (9) 

where we have defined C — >/~N /2, and the integra- 
tion variables as k = (ro/7r)Ak. By changing the de- 
tection angles (9,<f>), it is possible to measure the quan- 
tity J Ak for values of Ak that lie within a circle of ra- 
dius Rq — 2k-^ro/TT. The integral ([9|) has to be sampled 
with a set of values of Jqq within the region of inte- 
gration defined by the square k x e[— 1, 1], and k y e[— 1, 1]. 
Therefore, a basic requirement to obtain Jgp is to 
chose rofci > V2tt, so that the integration region is con- 
tained within the circle Ro. Note that this implies that 
kh / ki a tt > V%, where kiatt is wavelength of the standing 
wave lasers. An additional homodyne detection of the 
quadrature V\ — iek(aj CQ — Oka), allows one to obtain the 
quantity Jgp — jj^' m ^ — Jp " , ~ m - ) ; which combined 

with §§§ can be used to measure the operator J^ 1 '" 1 ^. 
Correlations in 3D can be obtained provided that both 
the direction of photon detection, and the direction of 
the incident laser are tuned to scan the whole momen- 
tum space. Following the same argument as before, the 
condition to get spatial information is r^k^ > ^/3tt. 

Once the operators Jp' m ' are obtained as the inverse 
Fourier transforms of field quadratures, it is straight- 
forward to calculate their correlations. These spatial 
correlations are useful to characterize some interesting 
states, like for instance cluster states \4>)ci- The later 
are defined by a set of eigenvalue equations for the op- 



erator & = ai ® Ienei?h y) oi, such that K^\(f>) cl = 
±|0)cZ @, 01 1 where neigh(j) specifies the sites of all 
atoms I that interact with an atom j. With this def- 
inition at hand, a cluster state in a 2D lattice can 

be characterized by checking that the spatial correla- 
tion (j(™i,mi) J (»i,m 2 ) J (- i 3,™3) J («4,m 4 ) J (n 5 ,m5) ) ^ equal 

to ±1 when the operators J z are next neighbors of 
J x . These spatial correlations can be obtained by 
making an inverse Fourier transform of the quantity 
(J Ak V Ak V Ak3 J Ak *J Ak5 ), that may be measured by 
considering a y polarized laser and correlations of dif- 
ferent homodyne detections of quadratures. Then, the 
inverse transform can be made following the basic re- 
lation |9]), and the analogous relation that exists for 
J^' m '. This scheme can be readily extended to mea- 
sure many operator averages which characterize mag- 
netic quantum phases, like for example, the valence bond 
strength, (Jj ■ J/), where j, I are nearest neighbors. 

One of the main limitations of the method is the fact 
that we need a good detection efficiency. The reason is 
that only a few photons are emitted at each solid angle 
(see Fig. (JT])), because the measurement time T has to 
be small enough to ensure that only a few atoms produce 
scattering and the measured state is preserved. Also, just 
as in TOF imaging, the reconstruction of atomic corre- 
lation functions, requires to perform each measurement 
M times, with M large so that the quantity 1/vM that 
characterizes the statistical error for the measurement is 
small. Nevertheless, one important difference with re- 
spect to TOF is that there are no limitations regarding 
the laser shot noise. In TOF experiments, the atomic 
noise has to exceed the shot noise of the prove laser. 
Here, the fluctuations of the prove laser can be easily 
distinguished and eliminated, since the laser polarization 
is different from the polarization of the emitted photons 
that we measure. Since our model is valid for measuring 
atoms in a Mott state, a further source of error would be 
the appearance of a superfluid region in the borders of the 
atomic cloud. This problem can be overcome by focusing 
the scattering laser to the center of the ensemble. 

In summary, we have shown how to detect spin cor- 
relations of an atom lattice within a Mott insulator 
state without switching off the potential. The detection 
scheme is based on the fact that spin correlations in the 
momentum space are proportional to correlations of the 
photons that are emitted in an off-resonant scattering 
process. Using different photon detection techniques al- 
lows to measure different types of spin correlations that 
are useful to characterize certain many body states, like 
magnetic phases. A complete sampling of these correla- 
tions in the momentum space can be used to obtain spa- 
tial correlations, which are useful to characterize some 
other phases like cluster states. 
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